General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA TECHNICAL 

MEMORANDUM 


NASA TM X-62,006 


* 

X 


s 

H 

< 

c/5 

< 

Z 


5-/im LASER RADIATION FROM A CARBON MONOXIDE GASDYNAMIC EXPANSION 


by Robert L. McKenzie 
Ames Research Center 


Moffett Field, Calif. 94035 
October 1970 


CN 

o 

<0 

2 

ae 

O 


M71- 11668 


(ACCESSION NUMBER) 

4% 


(PAGES) 

i _ 

£ (NASA CR OR TMX OR AD NUMBER) 



(CATEGORY) 



5 -ym LASER RADIATION FROM A CARBON MONOXIDE 


GASDYNAMIC EXPANSION 
Robert L. McKenzie 
Ames Research Center 


SUMMARY 


Laser power at 5 ym has been obtained from the gasdynamic expansion of 
admixtures of carbon monoxide. Mixtures of CO and N2, diluted in some cases 
with argon, were heated by reflected shock waves and expanded through a 
supersonic nozzle containing an optically resonant cavity. The resulting 
laser radiation, at wavelengths between U.78 and 5 *^ ym, originated from a 
multitude of partially inverted CO vibrational states. The spectral and laser 
power characteristics were consistent with radiative gain predictions from 
numerical solutions of the master rate equations for the vibrational relaxa- 
tion of anharmonic oscillators. The experimental ratio of multimode laser 
power, normalized by the total enthalpy flow rate or power supplied to the 
reservoir, was equal to or greater than that obtained from similar CO2 
gasdynamic lasers. 

^ INTRODUCTION 


Although the vibrational relaxation of diatomic molecules has been 
studied for many years , only recent analyses have brought to light the effects 
of some important molecular details that can significantly influence the 
relaxation process. For example, in the past, most theoretical studies 
assumed an oscillator model containing evenly spaced (harmonic) vibrational 
quantum states and used the corresponding Landau-Teller 1 transition probabili- 
ties in master equations describing their rate of vibrational relaxation. The 
results led to a relaxation process that proceeded from initial thermal equi- 
librium through a continuous series of Boltzmann population distributions, 
often described by a changing vibrational temperature. Analyses based on 
those assumptions provide useful results when the behavior of only the lowest, 
nearly harmonic, vibrational levels are of interest. For example, harmonic 
oscillator models have been applied to CO2 gasdynamic expansions in which 
population inversions between the lowest excited vibrational states of the CO2 
asymmetric and symmetric stretch modes were predicted. 2 * 3 Recently, Treanor 
et al. 4 showed the small uneven or anharmonic spacing of vibrational energy 
states occurring in most molecules to cause the upper level nonequilibrium 
population distributions to become non-Boltzmann, particularly in thermody- 
namic circumstances where the exchange of vibrational energy between oscilla- 
tors is a dominant energy transfer mechanism. In the case of rapidly 
expanding flows from a reservoir in thermal equilibrium, the relaxation of 
oscillators in each anharmonic vibrational state occurs at a rate that 


depends, in part, on the energy spacing and results in an overpopulation of 
the upper states. For extreme cases, a net inversion may occur. 

The overpopulation of upper states , occurring in nonequilibrium thermal 
expansions of anharmonic diatomic gases, is the physical principle on which 
the CO gasdynamic laser described in this paper depends. Other CO lasers have 
been reported previously in which nonequilibrium vibrational excitation is 
obtained from electric discharges 5-7 or chemical reactions . 8 The power avail- 
able from discharge or chemical lasers is limited by the size of the discharge 
or reaction zone and the rate at which molecules in the lower laser states can 
be removed or re-excited. Similar limitations are present in gasdynamic 
lasers, but they are less restrictive. Thermally induced inversions can be 
generated in large volumes by using large or cascaded nozzles and require only 
that the nozzle reservoir be heated to thermal equilibrium. Moreover, the 
expended species in the lower laser states are rapidly removed from the laser 
cavity by the supersonic flow. The available gasdynamic laser power then 
depends mainly on the mass-flow rate of heated gases that can be provided from 
any of a variety of energy sources. The power levels thus obtainable can 
exceed those possible in discharge lasers by several orders of magnitude. 

In the sections to follow, a theoretical analysis 9 of the vibrational 
relaxation and radiative gain in anharmonic oscillator expansions is first 
reviewed and provides a description of the physical details on which the oper- 
ation of diatomic gasdynamic lasers depends. The analytical results also 
served as a guide to the choice of proper experimental conditions and allowed 
a selection of extreme gas mixtures to be made. The configuration and instru- 
mentation of a shock-driven gasdynamic laser, similar to some CO2 devices 
reported previously , 10 * 11 are described next, followed by the experimental 
performance with mixtures of CO. Finally, a brief discussion is included 
describing some unique features of diatomic gasdynamic lasers and their 
comparison to the characteristics of similar polyatomic (CO2) systems. 


VIBRATIONAL RELAXATION OF ANHARMONIC OSCILLATORS IN EXPANDING FLOWS 


In their studies of the relaxation of anharmonic ally spaced vibrational 
quantum states, Treanor et al. 1 * wrote the master equations for the population 
rate of each state. They included the usual terms accounting for collisional 
exchange of vibrational energy with the translational energy of the incident 
particle (V-T collisions), plus additional terms accounting for the exchange 
of vibrational energy between colliding oscillators (V-V collisions). They 
showed that in circumstances where the latter V-V terms dominate, the solution 
does not require Boltzmann population distributions but, to zeroth order, has 
the more general form 

(1) 


= N° exp(aV - e'"' /kT ) 
J J J 


y 


Here, N is the number density of oscillators of species j in vibrational 

J y 

quantum state V, E is the energy of that state above the ground state, T 

J 

is the translational temperature, and a is an additional time-dependent 
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parameter. For the Enharmonic Morse potential, the oscillator energy of 
quantum state V above the ground state is described by 

E I = k6 j v t 1 * e j (v + 1)1 


( 2 ) 


where 0 is a characteristic vibrational constant of species j, and e is 
J ¥ 

a small anharmonic coefficient (e << 0 j. In vibrationally cooling nonequi- 

J J 

librium processes, a takes on growing positive values and, for small T, 
approaches the condition oc * 0 /kT. The product aV in equation (l) then 

J V 

nullifies the large linear terms in E and the resulting non- Boltzmann over- 

J 

populations of the upper states predicted by equation (l) are determined 
mainly by the magnitude of the small anharmonic coefficient, e , in 
equation (2). " 


Recent "exact" numerical solutions to the complete master equations have 
been obtained for supersonic expansions of CO in mixture with other monatomic 
and diatomic gases. 9 The formulation of the rate equations for gas mixtures 
is similar to an approach reported by Bray 12 and are given with their related 
V-T and V-V transition probability expressions in appendix A. The numerical 
solutions reported here differ from those of Bray, however, who used a first- 
moment approximation to obtain solutions for a constant temperature process in 
N 2 . These solutions contain no significant approximations and rely on a 
finite difference scheme for the numerical integration along an arbitrary 
one-dimensional flow axis. 


An example of the results is shown in figure 1 for a CO expansion where 
the anharmonic CO population distribution at expansion area ratio A/A# = 200 
is compared with a harmonic oscillator solution and a distribution based on 
equation (l) for the same conditions. The anharmonic upper levels are signif- 
icantly overpopulated when compared to the results of a harmonic oscillator 
model and, at larger area ratios, eventually become inverted. The distribu- 
tion of the lowest levels is well described by the zeroth order approximation 
(eq. (l)) as predicted by Treanor et al. 4 


The addition of a second vibrationally relaxing species, with energy 
levels nearly resonant to those of CO, can enhance the upper level populations 
through additional V-V energy transfer. In particular, N 2 is suitable because 
it has slower V-T relaxation rates than CO, allowing N 2 to serve as a reserve 
of vibrational energy (in expansions) which may be efficiently (rapidly) trans- 
ferred to the vibrational states of CO. Nitrogen therefore plays the same 
role with CO as with the asymmetric-stretch mode of C0 2 in all C0 2 -N 2 lasers. 
The vibrational exchange between CO and N 2 does not depend on the anharmonic- 
ity of either species, however, although it can be strongly affected by 
anharmonicity . Teare et al. 13 derive an expression for the locally equili- 


brated V-V exchange between harmonic CO 
temperatures , Ty : 


‘Vco 0 


N' 


and N 2 , in terms of vibrational 

e C0 

e N 2 " e C0 


(3) 


l v n . 
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For sufficient vibrational excitation (large 



translational temperatures, the CO vibrational temperature or energy can 
become large and, in some cases, exceed its initial or reservoir value. An 
example of this vibrational pumping effect has been demonstrated experimen- 
tally 13 and is illustrated in figure 2 by some numerical examples for anhar- 
monic CO and N 2 . The influence of oscillator anharmonicity is to invalidate 
the accuracy of equation (3) for large area ratios. Figure 3 indicates the 
range of validity of equation (3) for these examples by comparing ( T VC0^NUM* 

the numerically computed vibrational temperatures of enharmonic CO based on 


the relative population of the first and ground state levels, with (^Vco^Eq(3) 
computed with equation (3) using local values of T and Ty^ from the numeri- 


cal solution. The net effect of N 2 on the population distributions of anhar- 
monic CO is shown in figure 4, where a comparison is made with the previous 
100-percent CO case. The 5-percent CO — 95-percent N 2 mixture obtains a 
greater number of CO molecules in all levels above the fourth and displays a 
flatter distribution, both of which will contribute to greater radiative gain 
for transitions among those levels. 


From the preceding, one may conclude that for expanding flows the effects 
of anharmonic energy states on the vibrational relaxation process must be 
included if populations of any (except the lowest) vibrational levels are of 
interest. Also, the addition of a second vibrationally resonant species can 
significantly enhance the overpopulation of upper vibrational states in the 
primary species. The latter effect leads to a similar enhancement of radia- 
tive gain, a feature discussed in detail in the next section. 


INFRARED RADIATIVE GAIN IN DIATOMIC GAS EXPANSIONS 


The formulation of radiative gain coefficients for Doppler broadened 
vibration-rotation transitions has been summarized by Patel 5 in connection 
with a CO electric discharge laser analysis. The molecular circumstances in 
gasdynamic lasers are slightly different, however, because they usually oper- 
ate at higher pressures in the laser cavity. The combined effects of Doppler 
and collision broadening must then be considered. The formulation of radia- 
tive gain with both broadening mechanisms included is summarized in appendix 
B. The results are combined here with the previous vibrational relaxation 
calculations to provide an analysis of the radiative properties in diatomic 
gas dynami c expans i ons . 


Following Patel, 5 some qualitative features of the P-branch vibration- 
rotation gain coefficients with only Doppler broadening are shown in figure 5 
for a coefficient normalized by all of the multiplying molecular constants, 
which are independent of the upper level rotational quantum number, J, or 
rotational temperature, T r . Note that the vibrational population ratio 


xW 1 = 0.8, used in figure 5* 
J J 

sion; although as T r decreases 


does not represent a net vibrational inver- 
, positive gains appear and strengthen. In 
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the calculations to follow, gain coefficients are presented only for the 
P-branch rotational transition of maximum gain from each vibrational level. 

Predicted gain coefficients are shown in figure 6 for several expanding 
CO-N 2 mixtures. As expected, expansions containing N 2 yield considerably 
higher gains than those with CO alone. Since optical loss coefficients for a 
moderate size research device are typically near 10" 3 cm" 1 , the results of 
figure 6 indicate that at the conditions assumed, CO expansions without N 2 can 
be made to lase only in unusually low-loss cavities and afford no advantages 
over those containing large fractions of N 2 where much greater gain is 
achieved. 

A second aspect of the radiative properties from diatomic expansions is 
found in the gain enhancement due to reduced rotational temperatures as fig- 
ure 5 indicates. Since rotational relaxation requires only a small number of 
collisions, the rotational energy (temperature) of most expansions may be 
assumed in thermal equilibrium with the local translational temperature. 

Thus, to enhance gain, the local translational temperature at a given geomet- 
ric location in an expansion should be lowered without reducing the vibra- 
tional energy per molecule provided in the reservoir. This may be done 
through strictly thermodynamic effects by adding large fractions of a mon- 
atomic gas. Argon is particularly suitable for this purpose because, as a 
heavy collision partner, it also reduces the rate of vibrational energy loss 
to translation through V-T collisions. The enhancement of gain by diluting 
a 5-percent CO — 95-percent N 2 mixture with 80-percent Ar is shown in figure 7. 
At area ratios less than 300, an order-of-magnitude increase in the gain coef- 
ficient is predicted and maximum gain occurs at transitions from lower 
vibrational levels. 

The preceding results suggest three extreme gas mixtures to demonstrate 
the aspects of CO gasdynamic laser operation, namely, 100-percent CO, which 
is not expected to have sufficient gain except at very large area ratios 
(A/A # >> 200); 25-percent CO — 75-percent N 2 , which yields ample gain and has 
the possible advantage of having the greatest system efficiency because the 
number of additional species not directly contributing energy to the laser 
states is minimized; and 80-percent Ar — 5-percent CO — i5-percent N 2 , which 
should demonstrate the enhancement of gain over the previous mixture at lower 
area ratios and on radiative transitions from lower vibrational levels. These 
three gas mixtures were chosen for the experiments described in the following 
sections . 


EXPERIMENTAL ARRANGEMENT 


Mixtures of CO in N 2 and Ar were heated and supersonically expanded in 
the shock-tunnel arrangement illustrated in figure 8. The facility is oper- 
ated by filling the driven tube with a test gas to pressures near 1 atm. A 
shock wave is generated in the tube by rupturing diaphragms separating it from 
a high-pressure (200 atm) helium driver. The incident wave reflects from the 
driven-tube end wall, leaving a stagnant volume of compressed and heated test 
gas that acts as the nozzle reservoir and remains nearly undisturbed for 
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several milliseconds. The compressed reservoir gas ruptures a thin copper 
diaphragm and exhausts into the evacuated nozzle and dump tank. Expansion 
waves and helium from the driver then arrive and destroy the steady reservoir 
conditions. The nozzle was an axisymmetric 10° (half-angle) cone with a 
1.27 cm 2 sonic throat area, A # . Mirror ports, with their optical axis perpen- 
dicular to the flow axis, were located at geometrical area ratios, A/A # , of 
577, 1^55, and 2730. An optical cavity could be provided in any one of the 
ports by installing two circular mirrors, each with a 4.4-cm-diameter aperture 
and an inside spherical surface of 10-m radius. Power was removed through one 
mirror, made from a germanium substrate and coated for 3-percent transmission 
over the spectral region between 5 and 6 urn. Transmissivity rose through 
lU percent at U.7 ym. The other mirror was coated for maximum reflectivity 
over the same spectral region. 

Laser power was determined during the 2 ms period of quasi-steady flow by 
monitoring the temporal profile of the intensity pulse and measuring its total 
energy. 11 The intensity pulse was monitored with a fast-response Ge:Au detec- 
tor viewing an attenuated fraction of the beam reflected from a CaF 2 beam 
splitter. The transmitted laser power (93 percent) was absorbed by a cali- 
brated copper slug calorimeter. The combined pulse shape and energy measure- 
ments provided absolute time-dependent power values with an experimental 
accuracy of at least *20 percent at 100 W. The minimum power that could be 
resolved was less than 5 W. 

Spectral measurements were made with a wi de-band, grating monochrometer 
set at a different central wavelength for each run. The collimated laser 
beam illuminated a NaCl diffusing plate located ahead of the entrance slit to 
help fill the instrument optics. A 75-groove/mm grating, blazed at 8 yin, was 
used in second order with a 1-mm exit slit, giving a measured band pass of 
0.08 ym. The 5-vm shock-tunnel spectral signals were detected by a Ge:Au 
detector with a 1-mm-square cold aperture and a uniform spectral response 
between U.5 and 6 ym. The monochrometer band-pass measurements could not be 
made using the Ge:Au detector, however, and the shock-tunnel data showed evi- 
dence that a smaller band pass existed, apparently limited by the external 
detector optics. Hence, all spectral features may not appear in the experi- 
mental results although the major features not found initially were located 
by additional intermediate wavelength settings. 

Since hydrogen-bearing impurities are believed to have a large effect on 
vibrational relaxation rates, particularly in CO expansions, 14 attention was 
given to the test gas purity. Matheson "ultra-high” purity CO and Ar and 
"prepurified" N 2 were mixed, as received, in an evacuated and outgassed mixing 
tank several days prior to their use. The driven tube was also purged with N 2 
and evacuated before loading. These procedures are not sufficient to achieve 
the slowest vibrational relaxation rates reported by von Rosenberg et al., 14 
but represent a practical limit for convenient gasdynamic laser operation. 

From von Rosenberg’s data, 14 the conditions of these experiments are assumed 
to yield CO vibrational relaxation rates up to 100 times faster than the 
slowest possible and consequently have an adverse effect on radiative gain. 
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EXPERIMENTAL RESULTS 


Of the gas mixtures used (100-percent CO, 25-percent CO — 75-percent No, 
and 5-percent CO — 15-percent, N 2 — 80-percent Ar), only those containing N? 
resulted in laser operation. Examples of the instrumentation response to runs 
with the successful mixtures are shown in figure 9. The reservoir radiation 
(upper traces) was monitored with an InSb detector, through a 2-pm, wide-band, 
filter and quartz window located in the driven-tube side wall near the nozzle 
entrance. The reservoir radiation traces provide an indication of the com- 
bined variations of reservoir temperature and pressure. The reservoir pres- 
sures (middle traces) were measured with a calibrated piezoelectric transducer 
in a similar location. They indicate the arrival of the incident shock wave 
(small, step), followed by a large pressure step behind the reflected wave. 
Nearly-constant reservoir pressures were obtained for 2 to 3 ms in most cases, 
even for the off-tailored conditions of figure 9(b). The nozzle required 
about 0.5 ms to start as the delayed laser intensity in the lower traces indi- 
cate. Laser radiation originated from the starting shock wave and varied 
afterward in accordance with the reservoir conditions. The beam apparently 
contained a multitude of transverse modes since the beam size was between 1.5 
and 2.0 cm, 10 times the predicted single-mode size. Experiments with only 
the output mirror installed gave no signal in the Ge:Au monitor, verifying 
that the detected radiation came only from intense laser transitions. Note 
that while the laser radiation from the BO-percent Ar mixture (fig. 9(h)) 
decays soon after arrival of disturbances in the reservoir, the mixture with 
only CO and N 2 (fig. 9(a)) continues tc iase until most of the reservoir is 
exhausted. The extended period of the latter mixture is partly explained by a 
comparison of the spectral intensities in figure 10. 

The spectral intensities at all wavelengths from the Bo-percent Ar mix- 
tures were found to decay within 11 ms after starting (fig. 10(a)). The 
25-percent CO — 75-percent N 2 mixture, however, exhibited a first peak concur- 
rent with the multimode intensity, followed by a second more intense peak for 
some wavelengths. The second peak consistently occurred 11 ms after starting. 
The double spectral behavior of the 25-percent CO — 75-percent N 2 mixture is 
believed to be caused by the mixing of driver He with the test gas. Since lie 
is a monatomic gas, it also enhances the gain, although it increases the rate 
of vibrational energy transfer to translation. The adverse effects of Ho 
Icould not have teen too great, however, since the total multimode intensity at 
11 ms is reduced mainly in accordance with the reduced reservoir pressure. 

The spectral characteristics of all the mixtures at A/A# = 2f30, including 
those with the unknown He content, are summarized in figure 11. The C0-N2-He 
data represent the peak intensities 11 ms after starting, while the others 
were measured at 2 ms . In all cases, the intensities at each wavelength have 
been normalized by the maximum obtained at 2 ms and by the small variations in 
total multimode intensity from run to run. The spectral regions of eacn vibra^ 
tional transition are also indicated in figure 11 for the P(U) to P(l2) rota- 
tional lines. At A/A# = 2730, the 25-percent CO — 75-percent N 2 mixture lases 
most intensely near 5*05 pm, suggesting the transition of maximum gain to be 
from the sixth vibrational level. After 11 ms, the 5.05-pm radiation is com- 
pletely quenched and the shorter wavelengths are enhanced, in accordance with 
the effects of a monatomic additive (viz., He). The most intense radiation 
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from the 80-percent Ar mixture was near 4.78 pm, suggesting the transition of 
maximum gain to "be from the second vibrational level. Laser operation at 
lower wavelengths was suppressed, however, by the increased reflection losses 
of the mirror coatings below 5 pm (l4 percent at 4.7 pm) and possibly by 
strong boundary- layer absorption of the 1,0 vibrational transitions. The 
spectral data in figure 11 are corroborated by the gain predictions in fig- 
ure 12. Those calculations were done for the conditions of the experiment but 
used V-T rates 100 times faster than the established postshock wave values for 
both CO and N 2 . The calculations may therefore underpredict gain but are more 
compatible with the vibrational relaxation data for CO expansions reported by 
Russo 15 and von Rosenberg et al. 1 ** A comparison of figures 11 and 12 shows 
the spectral behavior of the laser to be consistent with the predicted transi- 
tions of maximum gain and gain distribution over the vibrational states. In 
all cases, the predicted gain corresponds to only partial inversions of the 
vibrational states and increases with area ratio. 


The total multimode laser power (from all wavelengths combined) varied 
with the expansion area ratio as figure 13 shows. The data reflect the pre~ 
dieted increase in gain with area ratio but are distorted by the effect of a 
decreasing fraction of flow passing between the laser mirrors as area ratio 
increases; that is, at the larger area ratios, most of the flow passes above 
and below the internal beam where its latent power cannot be extracted. The 
data may be compared with the projected performance of other nozzle configura- 
tions (e.g, , two-dimensional nozzles where most of the flow passes through the 
internal beam) and with the measured performance of some similar C0 2 gas dy- 
namic lasers 11 9 16 by recasting it. in terms of an efficiency based on the total 
enthalpy flow rate through the internal beam. This efficiency would then 
eqUol the ratio of laser power obtainable in a properly designed nozzle to the 
tot power required to heat the nozzle reservoir gases. For these 
experiments , efficiency is defined by 


Eff = 


P _A 
H 0 m LD 



where P is the measured jlaser power, H 0 is the reservoir specific enthalpy, 
m is the mass-flow rate through the nozzle, A is the nozzle cross section 
containing the optical axis , and L and D are the .Length and average diameter 
of the internal laser beam intersected by the flow. Efficiency percentages, 
computed from the power data in figure 13, are compared in figure l4 with the 
optimum output of two C0 2 gasdynarric lasers operated at Ames Research Center.* 
On this basis, the nonoptimized CO gasdynamic laser reported here is seen to 
match or exceed the capabilities of similar CO 2 devices. All of the lasers 
represented in figure l4 operated at reservoir enthalpy flow rates, H 0 m, 
between 1 and 3 MW. 

The influence of varying reservoir conditions is shown in figures 15 and 
l6. The maximum reservoir pressure (170 atm) was limited by the apparatus and 
no upper power level was achieved. The optimum reservoir temperature, for the 
80-percent Ar mixture (computed from incident shock speed measurements), 
appears to be near 2000° K. • 


Unpublished data taken from an arc-heated CO 2 gasdynamic laser by G. Lee. 



A second output mirror, coated for 10-percent transmission between 5 and 
6 wm, also supported laser operation at A/A # = lU 5 5 and gave approximately 
the same output power. The ideal mirror transmission will depend on the noz- 
zle configuration, area ratio, and desired spectral characteristics, however, 
and was not investigated further. 


SUMMARY AND APPLICATION CONSIDERATIONS 


The preceding experiments have provided evidence that supports the quali- 
tative nature of predictions based on recent anharmonic oscillator vibrational 
relaxation theories. Laser radiation, which depends on partial vibrational 
inversions induced in the thermal expansion of anharmonic CO, has been mea- 
sured at wavelengths corresponding to a multitude of lower vibrational trans- 
itions. Multimode power levels were obtained that match or exceed the 
performance of similar CO 2 gasdynamic lasers. 

Some unique features are provided by diatomic gasdynamic lasers operating 
on the vibrational energy states of a single vibrational mode. For example, 
laser transitions occurring between the upper vibrational levels, V and V-l, 
cause an enhanced inversion h »- tween the V-l and V-2 levels. Likewise, the 
relative populations between V+l and V are also further inverted. Thus, the 
removal of laser power from one location in a CO expansion should enhance the 
gain at several diffe*ent wavelengths farther downstream. This behavior is in 
contrast to the performance of a CO 2 system operating on the lowest excited 
levels of two vibrational modes. There, the maximum efficiency or laser power 
must be obtained by using cavity mirrors extended in the flow direction. The 
laser power is then rcntrolled by the collisional transfer rate of vibrational 
energy from N 2 to the CO 2 asymmetric-stretch mode (the upper laser level). 

The CO gasdynamic laser may therefore have the advantage of either providing 
equal amounts of laser power from laser cavities with shorter lengths in the 
flow direction or allow single-mode laser power for a multitude of wavelengths, 
with each wavelength from a cavity at a different location in a single expan- 
sion. These features accompany the disadvantages associated with the 
increased atmosphere absorption at 5 ym compared to that at 10.6 ym and the 
required disposition of large amounts of toxic CO from the nozzle exhaust. 


Ames Research Center 

National Aeronautics 
Moffett Field, 


and Space Administration 
Calif., 9^035, Oct. 1, 1970 
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APPENDIX A 


MASTER RATE EQUATIONS 


RATE EQUATIONS 


The number density of oscillators, N , of species J in vibrational 

J . ,.V 

quantum state V may be expressed in terms of a mole fraction, X , relative 
to the entire gas mixture, that is, 


X v „ 

J 


H 


1 


J.v J 


In general terms , the temporal rate of change of 
relaxation process may be written 


dX 

dt 


= I 2 ik« 

k 


V 

Jk^Jk 


(Al) 


due to a oollisionm 


( A2 ) 


where Q is the molar production per collision of X due to encounters 


\Jk 


J 


between species j and k, and Z is an average collision rate given by 

Jk 


Z = £ a 2 
Jk m Jk 



(A3) 


In equation (A3), P and m are the average fluid density and molecular weight, 

d 2 is a collision cross section for encounters between species J and k, nnd 

Jk 

u., is their reduced mass. If we denote the total mole fraction of species 
Jk ^ t 

k as = l X , including all vibrational quantum states L, and consider 
L k 

only single quantum transitions, the production term due to vibration-to- 
tranclation (V-T) energy exchanges becomes 




(A4a) 


The vibration-to- vibration (V-V) energy transfer contributes 
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r Q v l = P M ; 

r j? 

Ljkj v _v j 

[ ( Jk\v,V-l/ [ 


r -l L + 1 

I X k 


xp(E^ +1 ’ L /kT - E^’ V -‘/kT) - XjX^J 


+P 


jk 




x v+1 x" - X^' L+1 

J k J 


x£ +1 exp(E^ +1,L /kT - Ej + 1,V /kT)j| 


and the total production term is then 

4 ■ [4L * KL 


v-v 


(Ai+b) 


( A4c ) 


Equations (AU) include the detailed balance relations, exempxified by 




and are written in terms of relative transition probabilities 


p r * S 

= P r ’ 

jk 

Jk 

liM 

= P., 

\r,s/ 

Jk 


= P /p 5l ) 

jk\r ,s/ / jk\l ,0/ 


(A5) 


for later convenience. 

Terms similar to E^ ,V 1 in equations {Ah) are the energy difference 

between two adjacent states, V and V-l, of species j. They may be evaluated 
by assuming a Morse oscillator model in which the vibrational energy above the 
ground state is given by 


= ke^Vtl - e (V + 1)] 


(A6a) 


Here, 0. is a characteristic vibrational temperature and c is the small 

j V V-l ^ 

anharmonic coefficient. The energy difference, E * is then 

j 


P V,V-1 . 

J J " “~J 


E w . ,w ' = k© 4 (l - 2e .V ) 


(A6b) 
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TRANSITION PROBABILITIES 


Following Bray, 12 we can adopt a convenient form of the theory of 
Schwartz et al., 16 in which the probability ratios are given by 


— V , V- 1 
Jk 




(ATa) 


(ATb) 


The integral expressions of the adiabaticity factors F and F ’ are evalu- 

J J K 

ated by an empirical fit of Keck and Carrier 17 which bridges the gap between 
impulsive and adiabatic energy exchange, namely, 



e -2/3X) e -2/3X 


(A8) 


In equation (A8) , X is given by 


V |AE jk 



(A9) 


where Jij is a range parameter describing the intramolecular potential and 


AE 


Jk 


is the magnitude of the net vibrational energy difference for the 


exchange in question. Equations (AU) also require absolute values of the 
ground state probabilities, and Pj^i’q)* since > in many cases, existing 

theories do not successfully reproduce the experimental results, preference 

is given here to the experimental values in an empirical form. Data are 

available, at least for the pj»° values, obtained from measurements of the 

jk 

net relaxation time, behind the shock waves. Under the conditions of 

shock-wave experiments, Xj^ depends principally on P^° alone, allowing the 
latter to be determined with the correlation 


B., 


T Jk 


= exp(C Jk T- 1/3 ) 


(A10) 


and 


1.0 


ji 


_ T Jk p d jk *k^/w^kT 


(1 - e'V T ) 


) 


- 1 


(All) 
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Values of the empirical constants, B and C , used in this study were 

jK jK 

determined from the data of Millikan, White, and Hooker. 18 * 18 

Experimental data for the V-V probability P ^ ( i * q ^ are more sparse. 

However, what little is pertinent to this work (e.g. , refs. 20 and 21) shows 
that convenient closed- form expressions like those of Schwartz et al. 16 are 

usually suitable for predicting the temperature dependence of but 

JK A » u 

must often be adjusted in magnitude. We shall adopt that point of view here 

and obtain expressions for * q ) from the following: Schwartz et al. 16 

describe the transition probabilities for V-V transfer between energy levels - 
far from resonance - as 


P ( 8 * 1 ) 

p jk^.o' 


u 2 (AE ) 2 a 3/2 e"°J k 

o 39,2 (A V °J k 

^ M.M, „1„1 


E ?k 


1 - e 


- 2/30 


Jk 


(A12a) 


where 




r 

C 

= 3 - 


2 lT 4 y jk £ jk (AE jk) 2 l AE 


h 2 kT 




1/3 


\5 k 


2kT 


(A12b) 


and 


iE jk = lEj - 4 


In equations (A12), Mj is the reduced mass of the oscillator j, and 

Ej is its energy given by equation (A6) for V = 1. For the near-resonant 

case, Treanor 22 reports a modified version of the theory in reference 16 for 
small energy differences which reduces to the theory of reference l6 for exact 
resonance, namely. 


P (°» 1 ) - 

P jk^,0 J - 


p„.h 2 kT e _!2 J k 


Jk 




(A13a) 


where 


o 8 tt 14 0 2 

jk 25 jk J jk 


< aE Jk )2 

h 2 kT 


(A13b) 


As Treanor 22 points out, when either equations (A12) or (A13) are used away 
from the region in which they apply, the prediction is too small. Later com- 
parisons to the data of reference 21 showed values of 1/10 the larger of the 
results given by equations (B12) and (B13) to provide a reasonable estimate. 
Those were the values used here. The appropriate molecular constants are 
listed in table 1. 
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TABLE 1.- TRANSITION PROBABILITY CONSTANTS 


Species 

9j (°K) 

e j 

I? ( cm 2 ) 

(cm) 

CO 

3122 

6. 203x10” 3 

1.36x10“ 15 

2. 00x10" 9 

n 2 

3395 

6. 217x10“ 3 

1.38x10“ 15 

2.03*10“ 9 

Ar 

- 

— 

6.75 x 10” 16 

2.00xl0“ 9 


Ik 


APPENDIX B 


RADIATIVE GAIN COEFFICIENTS FOR THE VIBRATION-ROTATION 
BANDS OF DIATOMIC MOLECULES 


Radiative properties for laser applications are conveniently described in 
terms of a small-signal spectral gain coefficient defined by 


k(v) = - Un I v ) 


(Bl) 


where z is the optical path coordinate and I v is the spectral intensity at 
frequency, v. Since laser oscillation preferentially seeks the frequency of 
maximum gain (the line center) and gain measurements are most often made with 
the spontaneous emission nulled, the expression for gain coefficient in this 
application refers only to induced transitions at the line-center frequency, 
v Q . In terms of the local molecular properties, we can then write 


k(v 0 ) ^u^ v o^uH n u£ ^ v o ^ 


* N 1 

6 U fl, 

1 - — 

e*VI 


(B2) 


where N u and N^ are the number densities for the upper and lower radiative 
states, g u and g^ are their degeneracies, B ^ is an Einstein coefficient for 

induced emission, and n u j^(^ 0 ) is the corresponding normalized line-shape func- 
tion, evaluated at the line-center frequency, v . For purely Doppler broad- 
ened lines, usually considered in most laser applications, Penner 23 shows that 



where the Doppler line width is 


iv D 




tn 2 — 
m 


(B3) 


(B4 ) 


At most of the gasdynamic laser conditions, collisional broadening also 
becomes significant. The combined collision and Doppler broadened line pro^ 
files, evaluated at the line center, are shown by Penner 23 to be 


ex2[l - erf(x)] 


(B5) 


where x = -Win 2(Av c /A\>[> ) and the collision line width for species j is 


Av„ = 


- 1 £ 


c " it m f 


I *k«. 


AnkA' / 2 

V u jv 


(B6) 
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Values of the optical cross section, are listed in table 2 and reference 

23 for the gases of interest here. 

Equations (B2) through (B6) may be specialized to the vibration-rotation 
spectra of a three-dimensional rigid rotor, anharmonic oscillator for 
vibration-rotation transitions; u -* l = (V , J ) (V-ljJ 1 !) where V and J are 

the principal vibration and rotation quantum numbers. In this case, the 
Einstein coefficients, B for degenerate rotation-vibration states with the 

selection rule, AJ = *1, are correlated by 


B U* = b ‘° v (27Tt) 

where 

0, AJ = -1 (R-branch) 

n = 

1, AJ = +1 (P-branch) 


(BT) 


Since the P-branch gain coefficients are always greater, only AJ = +1 trans- 
itions will be considered further. The molecular constant, Biq> is determined 
from values of the commonly measured Einstein coefficients for spontaneous 
emission, A at wavelength v through the reDation 



(B8) 


The ratio 8uN./8o N u e< l ua tion (B2) is evaluated by assuming that the 


rotational states are populated according to a locally equilibrated Boltzmann 

distribution. On that basis, we denote the mole fraction of molecules in the 

V T 

combined rotation-vibration state, V and J, as X » and define 

V r V J V 

X = i X * , which is identical to the vibrational mole fraction, X , previ- 

j J 


V J 

ously computed. If AE ’ is the difference in energy between the Jth 
rotational state of level V and the ground rotational state of level V, 
then 


X v,j _ X v exp(-&E V,J /kT) 
6 J Qj.(T) 


(B9) 


where Q r (T) is a rotational partition function given by 
energy difference is 


V J 

AE » = k0 r J( J 


+ 1) 


1 - 6 r (v + |) 


Qj. = T/e r 


and the 


(BIO) 
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In the above, 0 r is a characteristic rotational temperature defined in 
Herzberg's spectroscopic notation 24 as 0 r = hcB e /k and 6 r = a e /B e . The 
ratio, ** or fundamental P-branch transitions is then equal to 


6 J X 


.v-i,j+i 


r V-l 


s j+i x 


V,J 


exp 


||( AE V * J - i.E V_1 ,J+I ;/kTj 


(fill) 


Similarly, the upper level number density in equation (B2) may be written 


p . . X V exp(-AE V *^/kT) 

N u = " (2J + 1) 2 


(B12) 


T/9, 


The formulation given by equations (B2) to (B12) is sufficient to compute 

. . v 

k(v ) for AV = 1, P-branch transitions, in terms of p, T, and X . The neces- 
— J 

sary radiation constants are listed in table 2 for CO. 


TABLE 2.- CO RADIATIVE CONSTANTS 


Bio (cm 2 /erg sec) 0 r (°K) 


(cm 2 ):k=C0 k = N 2 k = Ar 


*6.02xl0 5 


2.778 9.05x10“ 3 5. 3 y 10“ 1 5 5.3*10“ 15 1.7x10“ 15 


*Based on A(v, j:v-l, J-l) = l,lh sec -1 for v = 1; j = 7 at 


v Q = 2169 cm” 1 (ref. 2U). 
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Figure 1.- Population distributions at A/A* = 200 for harmonic and anharmonic 
CO expanding in a two-dimensional supersonic nozzle with a 15° divergence 
half-angle and a 0.127-cra throat height. Reservoir conditions were 
P 0 = 100 atm and T 0 = 2000° K. 



Figure 2.- Vibrational temperatures normalized by the reservoir temperature, 

T 0l in expansions of anharmonic CO-N 2 mixtures. The vibrational temperatures 
are determined from the relative populations of the first and ground vibra- 
tional energy state of each species. (Nozzle conditions are the same as for 
fig. 1.) 
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Figure 3.- Comparison of the CO vibrational temperature for a harmonic 

oscillator model calculated from equation (3) and the anharmonic CO vibra- 
tional temperature of the first and ground state populations from numerical 
solutions. (Nozzle conditions are the same as fig. 1.) 
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Figure 4.- The effect of N 2 on the population distributions of CO at A/A* = 200. 

(Nozzle conditions are the same as for fig. 1.) 
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Figure 5.- Normalized P-branch gain coefficients for Doppler-broadened 
vibration-rotation transitions from a partial {x v /x v ‘ ] < 1) vibrational 
inversion. 



Figure 6.- CO gain coefficients at A/A* = 200 in C0-N 2 expansions. (Nozzle 

conditions are the same as for fig. 1.) 
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Figure 7.- The effect of Ar dilution on CO gain coefficients in CO-N 2 -Ar 
expansions. V is the upper level vibrational quantum number. (Nozzle 
conditions are the same as in fig. 1.) 



Figure 8.- Shock tunnel and instrumentation arrangement. 
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(a) 25-percent CO— 75-percent N 2 , T 0 = 1800° K, P 0 = 100 atm. 
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(b) 5-percent CO— 15-percent N 2 — 80-percent Ar, T 0 = 2100° K, P 0 = 100 atm. 
Figure 9.- Reservoir and laser intensity instrumentation records at A/A* = 2730. 
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(b) X 0 = 5.14 ym 

Figure 10.- Spectral and total intensity profiles from 0.25 CO - 0.75 N 2 at 
A/A* = 2730, P 0 = 100 atm, T 0 = 1800° K. Intensity level response of the 
two detectors is not the same. 
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Figure 11.- Spectral characteristics at A/ A* = 2730. Lines connecting the data 

points may not indicate all spectral features. 



Figure 12.- Gain predictions for the experimental conditions. 
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Figure 13.- Multimode power at 2 ms after starting. The boundary- layer 

corrections apply to the 25-percent CO - 75-percent N 2 data only. They are 
the result of some unpublished measurements made in N 2 flows at the same 
nozzle conditions by G. P. Menees. 
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Figure 14.- Power efficiency based on equation (4). Boundary- layer corrections 

are not included. 
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Figure 15.- The effect of reservoir pressure on multimode power at A/ A* = 1455 

and 2730, T 0 = 1800-2100° K. 
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Figure 16.- The effect of reservoir temperature on multimode power for 
0.05 CO - 0.15 N 2 - 0.80 Ar at A/ A* = 2730, P 0 = 100 atm. 
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